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FRACTIONAL RELEASE OF A TRACER ELEMENT 
THROUGH A MOVING BOUNDARY 

Stephen D. Lowe' 1 ' 

Lawrence Radiation Laboratory 
University of California 
Be rkeley, California 

October 30, 1963 
ABSTRACT 



A full solution for the fractional release of a tracer element 
through a moving boundary has been obtained for the cases of a slab 
and a sphere. In both cases, the initial distribution of the tracer is 
assumed to be uniform throughout the body. For a constant rate of 
boundary motion, the full solution is valid for all times over the range 
from zero to complete evaporation of the body, and this solution may 
be applied to all finite rates of boundary motion. 

By applying certain limitations to the time range and to the total 
amount of boundary motion, the full solutions are reduced to approxi- 
mate forms that are presented graphically as a family of curves. 

For a more limited range of application, i. e. , for small values 
of the quantity (~t) (< 0.2), where b is the rate of boundary 

motion, t is the time, D is the diffusion constant, and a is a charac- 
teristic dimension of the body (half-width for the slab and radius for 
the sphere), simplified expressions are developed for the fractional 
release of the tracer, f, which may be used to more accurately deter- 
mine the diffusion coefficient, D. 

These simplified expressions are 

2 (“) 1/2 + Kt ) for ,he slab - 



f = 



and 



f = 



IT 



7T 



1/2 



1/2 V 



^ty/2 + 3 f or the sphere. 



Lieutenant, Civil Engineer Corps, United States Navy. 
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I. INTRODUCTION 

A. Basis for the Problem 

During his investigation of the diffusion of xenon in uranium 
monocarbide, Shaked noted the slope of his plot of the fractional re - 

l/Z 

lease, f vs (t) ' increased noticeably during a high -tempe rature 

1 

anneal. He postulated that evaporation of the uranium monocarbide 
was responsible for the increase. 

l/2 

A similar increase in the slope of the f vs (t) 7 plot had been 
reported previously for experiments in the temperature range, 1600 
to 2200 a C. 2 

With the increasing interest in high -tempe rature applications 
for uranium fuels, in thermionics, for example, it is considered of 
interest to investigate the release of a fission product from the fuel 
element due to the additional process of evaporation occurring simul- 
taneously with the process of diffusion. 

B. Discussion of the Problem 

Several investigators have reported diffusion constants for 

3 

xenon- 133 and krypton-85 in uranium dioxide. A lesser number 

have reported results for the diffusion constant of xenon- 133 in 

\ 

uranium monocarbide. In the "conventional ft analysis, one assumes 
a uniform initial concentration of the tracer element in the geometri- 
cal body whose boundaries are fixed with time. Evaporation of the 
body is discussed in a general manner, but the effects are minimized 
experimentally.^ ^ 



C. Statement of the Problem 

This report is an investigation of the effect of evaporation upon 
the fractional release of a radioactive tracer element from a fuel 
body. As such, the decay of the tracer is neglected for the purposes 
of this report (i. e. , K = 0), although this effect can and should be 
taken into account for experiments involving long postirradiation 
anneal time. 

An expression for the combined fractional release due to both 
diffusion and evaporation (i. e. , diffusion through a moving boundary) 
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is obtained. This result is compared to the fractional release derived 
for the case of M pure diffusion 1 ' (i. e. , diffusion through a stationary- 
boundary). 

The analysis is performed for the case of the finite slab and the 
case of the finite sphere. It is assumed that the initial concentration 
of the tracer element in the body is uniform and that the concentration 
of the tracer element at the boundary (moving or stationary) of the 
body vanishes. 



II THE SLAB PROBLEM 



A. Concentration of the Tracer 



One seeks a solution for the diffusion equation 




_(a-tt)<x<a-£t (ii -1) 



such that 



C(x J o) = c o 



( II - la) 



and 



C L±(d-bt) l t~! = o -For o <t <°/b J 



(II- lb) 



where b is the rate of movement of the boundary that may be obtained 

from kinetic theory if one assumes evaporation in a vacuum of a 

4 

material of known vapor pressure. The assumption of a constant 
evaporation rate ignores the possible effect upon vapor pressure, and 
hence upon the evaporation rate, if there is a change of chemical com- 
position of the solid surface accompanying evaporation. At time 

t = a/b, the sample will have completely sublimed. 

5 

Chambre has shown that the solution to Eq 0 (II-l) that satisfies 
an initial even function f(£) and Eq. (Il-lb) may be expressed as 




(H-2) 




where for the case under study 




(Il-Za) 




(H-3) 
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Making the substitutions into Eq. (II-2), one obtains the result- 
ing expression for C(x, t) : 




(II -4) 



Cfat) Z(jrDt)' /s - 

m-i -o~ L -* 

Performing the integration over £ (see Appendix D), one may show 

+ f | f 0”ex P f K * - #)]} {exp {^)[erf^^ 

-erf 



— - 2U*)* -j; (II - 5) 

+ 2 (- if exp {- ~l££)jj {e X p(-iii-)|er-f 

erf 52(DtV* J j. 



This may be rewritten as 



■] 



Ill & II Id y UC J. C W J- J- x± 

/j /* Dt') - Co fC r -jT 1+Va. . er J? /- Vcl 

c (a > ^ ~ ^ L er+ ^m^ +en » , r , wv , 



J 



(II-6) 



+ 



zXWT* j j 

^ 2f 0 exp {exp [- s^ft)] 

^ l "a?/ 



Making the definitions 



_ (Zh-l) + (ya.-Zh^4jr) ll 

^r¥F^ J j • 



r- 



- vt 

c& 



(II-7) 
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and 



whe re 



n - 

r T> ’ 

0<f$Z 



(II-8) 



and substituting into Eq„ (II-6), one obtains 

mfr % fa 



(II -9) 



+ xp|- r,p [n(i-pq-^.]j jerf X/a - 

-erf^dd^d±2ik.| 

+ j-ltt-i)' exp [erf ^ ( £p)+ h * L 

_ er r Z„ 0 - 6 Z) r I +*/<,. ! 

n* j 

By using Eq, (C-4) and setting B = 0 (i. e„ , no evaporation or boundary 
motion), Eq. ( II -9 ) becomes 

c{\z) = c„- ^ [erfe(^)+erfof-d^) 



- & l(-/f 

/C nz, 



<rfc (Zn-fQ- X/a. _ 

Tz 7 ^ nc zz' /x - 



(II- 10) 



+ erf c &^-erf c & 4 £k]. 



This may be rewritten, after combining terms, as 



C(V) = C. - CZ(-iT [erfc ^)-^ + eHc- ^^ ] • 



Equation (II- 11) may be recognized as one equivalent to Eq, (A-6), 
which shows that the expression for the tracer c oncentration does re- 
duce to that for the stationary boundary case when no evaporation is 
conside red. 
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B. Total Fractional Release 

The total release of the tracer element due to the combined ef- 
fects of diffusion and evaporation may be obtained by integrating the 
concentration at any time (t) across the volume of the slab, dividing 
the result by the total initial amount of tracer present, and subtract- 
ing this quotient from 1. Thus, 



r = Q.-Qft) 

+ Qo 



I _ Q (?) 

' <So ' 



( 11 - 12 ) 



where 




(II- 12a) 




refore, n-^) 

_r _ i _ z!o c(y*>z)d(v*) 
T ~ 2 C 0 




- e 
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= i -*rv 



°"^ r - ^ + er f J- Vc 

(I-/3Z) 



ZZ 7 ^ 



zz 



' fe ] d (%.) 



(11-13) 



- z^f''^ ex p[-T^ ( '^}/«p[^j] (erf ^'W^-Vp. 



zz '^ 



- er p 2»O-0$-l- % 

- i i ( -f e *P W lerf M^tl^L 




_ er fMt0d^kU(^ . 



Performing the integration and combining, one has 

~f- f3Z + Z' lz [ierfc ~^r z - Cer-fc ] 



+ t'^ 



E n / t < n -'X , T*$ n **] ""'fftrf %y ' - 



-l”erfc &dXt£tl±L 
7.X'' 7 - 



+ 



T /z £ (-lT +W) exp f- V(3 (yn-$-p z)] finsz'^T ' 

fytozl [J o n-'/z. 



ZZ /z 



(11-14) 

ZZ^ 



For r — 0 ( t — 0) , the fractional release, given by Eq. (11-14), is 
zero, as is to be expected; and for (3t = 1 (or complete evaporation), 
the expression becomes 

+ 7T ^ 2t Vx ~~ ^ ^z" 7 ^ ] 

-/ 

' WTC 



-f 



+ i n 'erfc ^Jph. - i n erfc j , 



c 2t* " L '" er ^ £•*. 



which is equal to 1, after cancellation of equivalent terms, and cor- 
responds to the total release of the tracer. 

Setting p = 0 (no evaporation), Eq. (11-14) becomes 



■f 



= =4* t"*-r t Wc4’?, 



TT'/Z 
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(H-15) 



This may be reduced to 



~f=ZZ /Z [_ n ZSC' 1 )” J> 



( 11 - 16 ) 



which may be recognized as equivalent to Eq. (A-8), the fractional re- 
lease under conditions of a stationary boundary. 

C. Approximate Expression for the Fractional Release 

For the usual values of r encountered in fission gas -release 
experiments (t < 0.01), the expression obtained for the fractional re - 



If the further limitation is made that (3t< 0.5, Eq. (11-17) may be 
reduced further to 



lease in Eq. (11-14) may be immediately reduced to 







(11-18) 
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of which the n= 1 term is of the greatest significance, as follows: 

f- fi'Z + Z' /z terfc 

6 , (11-19) 

By using the relations (C-13) through (C-18), Eq. (11-19) be- 
come s 

-P=pz+ -‘'p'rerfc sg!' j 



- j3 f Z + ZZ' lz LCirfc '@=r~ l + 




( 11 - 20 ) 



. 1/2 



. 1/2 



For large (i. e. , >2.5), Eq. (11-20) reduces to f = |3 t, 

which indicates that the release of the tracer is due solely to evapora- 
tion. 

For p - 0, Eq. (11-20) becomes 



-f- 3±. 'r ■> 

■ ~ ItVtl <~ 



(II-21) 



the approximate expression for fractional release from a slab with 
stationary boundaries. 

Figures II. 1 and II. 2 show the fractional release from a slab vs 
[time] i/fZ (t 1 / 2 ) fo r various rates of evaporation (p). As a reference, 
the stationary -boundary plot is given, and in each figure it is the 
lower line. Figure II. 1 is for the longer periods of time, whereas 
Fig. II. 2 represents shorter periods of time. 

Figures II. 3 and II. 4 demonstrate the contribution of the individ- 
ual terms that make up the expression [Eq. (11-20)] for f. Figure 
II. 3 represents an evaporation rate of p = 10 for times less than 
r = 0.0l, and Fig. II. 4 covers the time spanless than j = 0.0001 for a 
higher evaporation rate of p = 50. 



-10 



(1) 


0 = 0 


(2) 


0 = 2 


(3) 


0 = 10 


(4) 


0 = 20 


(5) 


0 = 50 


16) 


0 = 100 


(7) 


0 = 200 



(no evaporation) 




MU-13935 



Fig. 



II. 1. Fractional release vs 
of evaporation. 



T 



1/2 



for various rates 
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(1) 0 = 0 (no evaporation) 

( 2 ) 0 = 20 

( 3 ) 0 = 100 

( 4 ) 0 = 200 

( 5 ) 0 = 500 

( 6 ) 0 = 1000 

( 7 ) 0 = 2000 




MU. 33926 



Fig. 



II. 2. Fractional release vs 
of evaporation. 



T 



1/2 



for various rates 
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f = pr+ 2 t 



1/2 



ierfc 



£I 1/2 + t 1 / 2 g (4) m ( , £!: l/2 ') m i m+1 erfc£l 

2 m=l \ 2 / 2 

— Release — no evaporation ((3 = 0) 

- - Combined curve for (3=10 



(1) Pr 

(2) ZT^ierfc^’ 72 

(3) T 1 / 2 (4)Ql 1/2 ^i 2 erfc 



Pt 1 / 2 
2 



(4) T 1 / 2 (4) 2 (^ l/2 Ji 3 er£c^ ,/2 

(5) T l/2 (4) 3^‘/2j.4 erfc ^l/2 




MU- 32927 



Fig. II. 3. Components of fractional -release curve. 
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pT + 2T>/ 2 ierfc^ V2 + T l/2 £ <4) m (?i l/2 ') m i m+ 1 erfc 
1 m= 1 \ z ) 

Release — no evaporation 

Combined curve for 0 = 50 



0t V z 
2 



( 1) |3r 

(2) 2 T l72 ierfc El ^ 

2 



(3) T '/ 2 ( 4)0I l/2 )i 2 erfcf l/2 

(4) T‘/ 2 (4) 2 (^ l/2 ) 2 i 3 er£cf: l/2 




Fig. II. 4. 



Components of fractional-release curve. 
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III. THE SPHERE PROBLEM 
A. Concentration of the Tracer 
A solution for the diffusion equation is sought. Let 



iC(r,t) 

TT - -J-’l ir' — 



Z ?C(f,t) J 
r ~ 7 ? J 



y 



0 <r<a-kt 



such that 



C(r,o)= C. 



and 



C 0 -for 



0<+ 



(III- 1) 
(III - la) 
(Ill-lb) 



where b, again, is the rate of boundary movement, At time t = a/b, 
the sample will have been completely sublimed. 

Equation ( III - 1 ) may be reduced to that for the slab problem by 
making the substitution 

u.(r } t) - r C( (ni-2) 



subject to the conditions 

u(r,o) = rC a 

and ,'k) ~ 0 -for 0<_' 



(III-2a) 
(III -2b) 



Since C(r, t) is finite at interior points, it is necessary to define 

u( o,-t) = o • ( ni - 2c > 

Chambre has shown the solution of Eq. ( III — 1), which satisfies an 
initial even function f(£) and Eq. (Ill- lb), to be 

CW) - { Ij^ ex P t Taf] 



( HI -3) 
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where, for this problem 

-ftj) = 



(III -3a) 



and s(r,tj$) = e*p [-%$-] 



" A r[ ~VDF 



(III -4) 



By making the substitutions into Eq. ( III - 3 ) , the resulting expression 
for C(r, t) is 

c (r ^ = y* fiJ 

A 



-6^ 



r , n (- (r-i-2rutf 

eK P[ 



(III-5) 



+ 2j$ex P [^(^] 

-„ p (-Si^£)]lP v 

By performing the integration over £ (see Appendix E), Eq. ( III - 5 ) 
may be shown to be 

c(M0= c„- ^ [«^^- erfc ^ 



_ C.a. 

X'r 



afoty'- ^ w*- J 



r 2n (a. -kb) - (a +t~) ~] 

+ er + c zmv^ — J 



+ 3^is exp[- g(n*.-M +!•)] [< srfc 



( III -6) 
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Using the definitions ( II - 7 ) and (II-8), then substituting into ( III - 6) , 
one obtains , . 

__ Co l-r/z. „._r. ! + r /a 



CM ' C €r ^ c ] 






l- r /c 



%. - te 




+ 



C 



ZZ ,/z ZZ' 

~ j$6) 1 & p l" 7 / 5 ^ Z) ' { erfc 

+ erfc ZnO-fld-l-Vo- 

zxP 1 

+ !+ 7cl 



X r /*) my 




£exp)-n 



ZZ J/ - 



fOu 



+ e^fc 2n(l-{3z)-l + r/ c 

il \ 

- fy~^ ex pf'Tt”('t T )' 7 d} 1 Cf> *£ l ~ l ' / °- 

- />^ O/aO-/- I i& - 

P T-'/t 

/ i ^ ^ 

+ c r /l 



(III -7) 



7a S 



exp j - ni?,[ r i(i-fiV+ r h}{ I Hrfc ’ l ' + r l'~ 

-Url 

Setting (3 = 0 (i. e. , no evaporation or boundary movement) in 
Eq. ( III - 7 ) , it follows that 

CM - C. 

-fe?| 

- Pt4r ^'h r la. _ a yX, &n-l) + r L 1 

2 Z' 1 * C 2 Z"i J 

- ^ r "* (W c - cede -^] 
-<firk &gfc$. y- i€rfc ( Z ^^ ] 



( HI ~ 8) 
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The ierfc terms reduce to zero, and after some recombination, 
Eq. ( III - 8 ) becomes 



This may be recognized as equivalent to Eq. (B-6), which shows the 
expression for the tracer concentration to reduce to that for the 
stationary -boundary case when no evaporation is considered. 



The total release of the tracer element due to the combined ef- 
fects of diffusion and evaporation may be obtained by integrating the 
concentration at any time (r) throughout the volume of the sphere, 
dividing this result by the total initial amount of tracer present, and 
subtracting this quotient from 1. Hence, 




B. Total Fractional Release 




(III-10) 



where 




(III-10a) 



and 




( III - 1 Ob ) 



Thus, 
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Integ rating and combining the terms in (III -11), one finds 

■f * 3(5 T - 3(^rf i-SZ'^Cl-pZ^Lerfc-^Sier-fc. 

~CZ (ferfc -^~, r L*trfc ) + 6 ZO-^Cerfc |£, A -f- 1 2 er/c 2 ^ 



-IZ X V * ( der/e If,. - t 3 erfc 

+ 3TZ'-^)ij ; in^T'e>< P + - ' 

+ rerfc &dXti^Ll 1 

^ £T //3 ~ J 

+ 3 Z lh (i-($z)Z(-0 nf \2fz3&p 

■+ rerTc (^t 'KlrAlbl 1 
^ ^7u L 

HZ 2(-<r ? 7 fex p [-hp(n-/Xi-pt)Jj^ r f c ^zl^±l 

h1=? 2.T ** 






■ C Z£(- ‘' 1 ') r 

y/)z£ /2 7T ( HI -12) 



+ (derfc Zd 



~ C J 

+ a (/ TK^f - '' (V "f ex p[-Y"-M] [iTtrfc. i^Ebl 

yy)^ 7 ^77 

- L^erfc &Z^SlH iJ 

HZ O-ptyZy >T( 2r f' 1 ' 1 ') e*P HfZMp-pV] lperfr^Z^hl31 

yy] - 2j /Co 

1 2 7 !/ 'J5 <-lj"(r^?Hf J % 'T exp [- t ’p( n it'-ftf] y_r-fc @”- l X l p£‘d f 1 

vtj=3 ^ 

- iVc^/j 

■l2T*2(-» m (^WyvTtxp[-n^^^ I 

yy> -j XL 23 z 



+ 
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For r = 0(t = 0), the fractional release, as given by Eq. (Ill- 12), is 
zero, as expected; and for |3 t = 1 (or complete evaporation), the ex- 
pression becomes 

f = 3-3*1 -CT(c^rfc ^-4- 

- IZ r 3/2 ( 3erfc 2 -Pfe -3erfc 
+ CZS(-r'(^f' 1 (^erfc^ * rerfc 



-/ 
YH-Z 



m - 3 



- rerfc ^ -rerfc 



(iW c ^ ^ 

>^=3 . / - ! V 

— i^er-fc 2% r h. > 



zz^r 

which is equal to 1, after cancellation of equivalent terms, and cor- 
responds to total release of the tracer. 

Setting (3=0 (no evaporation), Eq. ( III -12) becomes 

■f = ^ Z /z +■ 3 Z‘ /z (<erfc -f 7 ^~ -3 T -hZZ^erfc + v?- T 

+ CZi^rfc^- lzr Vi (i TT-'^tlZZ^ rerfc ^ 

+3t*| Wc 



± LZZ^ | i z erfc^^ + i z erfc 

- 6 ^er/y: i *crfc p j 

^-ilerfc^ + L'erf c ^ - i'erfc fijj 

-l22f /z £ 1 3-rk ^-Serfz&h - ferfc^ +l 3 erfc ^ j • 



( III -13) 



This may be reduced to 

- f = [ 7T "' /2+ 2£, Ceir fc &1 ~ 3Z, 



(III - 14) 
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which may be recognized as equivalent to Eq. (B-8), the fractional re- 
lease under conditions of a stationary boundary. 



C. Approximate Expression for the Fractional Release 

For the usual values of t encountered in fission-gas -release 
experiments (r < 0.01), the expression obtained for the fractional re- 
lease in Eq. (III-12) may be immediately reduced to 

+3Z /lz 0-f3z)ter-fc 

+ CZ(l-(3Z)l z crfcj$ lx - lZZ S/z C 3 e^c-ff 



T Z 0-fz)E( if j? X p [-np(n- ifif-z)] t z rfc 

p[-hp(n+iYi 1 *z)yv er f c fa+ifagzh l | 



+ ex 

t 



(III-15) 



+ f x p[- ridh-iYh-ar)] i z l 

vn-Z l 1 J 

- €X p [- rfc 1 

-6T(/-pz)Zf-/r(ln f }Z' , f' ! ppl-upfr-tfi-pz)! 

m=z ~2X' Zr~ 

expLM r ' H M'l,^ P^IXl-fz)- 



+ 



m-3 Zz ,'u “ 

— exp [-hftwYf-pr)] o^erfc ^dZJppllz ! I 

^ Z* 20 / 



By making the further limitation that (3 t < 0.5, one obtains 



f = 3fz-3((izf+<i3iT-i-3Z' ,z (l-{i'^^i'Xc jgr, 
- C.Z L’trfc^ + c Z0-{3z)t*er-& ^ 
-l2Z ZM c 3 erfc ^ 
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^ (m- i 6) 

+CZZ(-irt” , -ifc n p exp [-vftHXi-pzj] L^erfc. & n ~ 

fl-l < L 

Wl-2 




of which the n = 1 term is of the greatest significance, as follows: 

■f=3p?- 3(f3zf+(i szf + 3Z‘ J '(i- 

-iU*erfc^ * 6, Z(l-^) terfc ^ 

- I2Z 3/i l s erfc 
+ 3 Z 

m ~-' (III- 

+ 6T Ji' , ^- , ^/ 3T 'T'VeH^c 

- CZ(l-pz)g(-0*fepz'^) m ' Z C m erfc. 

- IZZ^i-irfa-zXZfs^y"* i”erfc ■ 

M-3 (C L 

By using the relations (C-13) through (C-19), Eq. ( III -17) becomes 

f = 3pt(t-(rc)+(pxf- 3T - 3Z0-(3z) 

+ & Z /z (i-p' t ) ieAc 4 — 4 + IZT(i- (3 t) i. er-fc -4^ h 
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+ 3 z k Ml w ’ tf-T * I 1 " 

(111-18) 

-tcTO-pt)l arypT i-%wc -£■" 

+ CZZ ^ Wr^Y T%rtc -££-* 

+ /2 r 3/ ^| i^T^'friAc ^ ■ 



The last three summations will be small compared to the first terms; 
therefore, f will be approximated by 

^3p'l(l-p‘i)+(pzT-3X-3Z(l-(3't) +IZZ(l-pz) 

( III - 19) 

+ 3't k 0-(x)farfcjZ VJ (VJ7- f^T'rkrfc ■ 



Bt 1/2 Bt 1/2 

For large (i.e., >2.5), Eq. ( III - 19) 

-f= 5fi'r(l-f3z)+(p‘Z?-3Z-3'Z(l-(3'z) 1 



will reduce to 



which indicates that the release of the tracer is due primarily to 
evaporation. 

For (3 = 0, Eq 0 ( III - 19 ) becomes 



C — C ? 

T Jf '/ z . C 




(III -20) 



the approximate expression for fractional release from a sphere whose 
boundaries are stationary. 

1/2 

If a further restriction is made, such that r <0.01 (or 
t < 0.0001), then those terms in Eq. ( III - 19 ) having t as a coefficient 
(excepting (3t terms) may be neglected, and one obtains 
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-f = 5 fSZ 0-(3l) +(jet) 

+ 3Z' h (hf3z) \ 2 tei-fr -^-'V g (hT^J C m kr-fc-gf k j • 



(III -21) 



«t 1/2 

For large 1- , Eq. (III-21) reduces to 



-f= 3p'C(l~(3'Z) + (/S'cf 



which indicates the release of the tracer to be due solely to evapora- 
tion. 

For (3 = 0, Eq. (III-21) becomes 




(III-22) 



the more familiar expression for fractional release from a sphere 
with stationary boundaries. 

Figures III. 1 and III. 2 display the fractional release from a 
sphere vs [time]^^ (r^^) for various rates of evaporation ((3). As 
a reference, the release for the case of a stationary boundary ((3 = 0) 
is given for each figure; see line (l). Figure III. 1 (note the dotted 
line) in addition, shows the effect of neglecting the 3 t term in 
Eq. ( III - 2 0 ) for longer periods of time. Figure III. 2 represents the 
release of the tracer element for shorter time periods. 



s 
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f = 




T 



1/2 



(1) 



(3 = 0 (no evaporation) 

1 = — — r — - 3 t 

TT 



(2) (3 = 2 



(3) 


p = 10 


(4) 


O 

<NJ 

II 

CQ. 


(5) 


p = 50 


(6) 


p = 100 




1/2 



Fig III. 1. Fractional release vs t 
of evaporation. 



for various rates 
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(1) 


p = 0 (no evaporation) 


(5) 


p = 500 


(2) 


p = 20 


(6) 


p = 1000 


(3) 


p = 100 


(7) 


p = 2000 


(4) 


P = 200 








MU- 32930 



1/2 



Fig. III. 2. Fractional release vs r 
of evaporation. 



for various rates 
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IV. DISCUSSION 

A. General Results 

For all values of (3>0, the fractional release with time is in- 
creased over the corresponding value for f when (3 = 0. For any value 
of time (or t), f will increase, for both the slab and sphere cases, 
as (3 increases. With smaller values of 6 , say (3< Z0, the curve of 

f vs very closely approximates a straight line for t^^< 0 . 01 , 

a/j 

which indicates that f is proportional (or very nearly so) to t 

l/z 

The same statement may be made for r <0.1, provided ( 3 < Z. 

B. Initial Release 

This theory does not account for the large early release rate in 

1 6 

postanneal experiments noted by several observers. ’ In these 

l/z 

cases, the slope of the f vs t ' curve decreased with increasing 
time until a region of more or less constant slope was attained. 

C. Spherical Case --Stationary Boundary 

Figure III. 1 illustrates the effect upon f of neglecting the second 

a /Z 

term of Eq. (III-ZO). For values of r < 0.01, the difference in 
values of f resulting from lack of consideration of the second term of 

i/z 

Eq. (III-ZO) is negligible. For larger values of r ' , the fractional 

release predicted by Eq. (Ill— ZZ) becomes increasingly greater than 
that predicted by Eq. (Ill -Z0), and differ s by 0.03 ( ~ 10%) at t^^ = 0.1. 
Since both Eq. (III-ZO) and Eq. (III-ZZ) appear to be plotted as straight 
lines, an error of 17% in the value of D results for the above value of 
r from neglecting this term. It is noted that Eq. (III-ZZ) closely 
represents the release curve resulting from boundary motion for 
which (3 — Z . 

D. Spherical Case --Moving Boundary 

i/z 

For values of t ' < 0.01, Eq. (III-Z1) may be used in preference 

to Eq. (Ill -19) for all values of (3 t . If (3r (remember p r independent 

of D and is bounded between 0 and 1 ) is approx<0.Z, then the term 
3 

(pr) may also be neglected. 
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1 l/z 

For large values of ^ p t 1 , the fractional release is approxi- 

mated by 3(3 t( 1 - (3t) +(|3t) 3, which is a cubic equation in time (or r). 
This may occur at short times for very rapid boundary motion (6), or 
for very long periods of time (t). Neither of these cases, however, 

is useful in making a determination of D. 

1 l/2 

For small values of — fir ' , pertaining both to relatively short 

periods of time and small rates of boundary motion, the fractional re - 

l/2 

lease is directly proportional (or very nearly so) to r 1 and there- 

1/2 

fore to t / . It should be possible, therefore, to develop a simple ex- 

pression that may be used to determine D with greater accuracy than 

would the use of Eq. ( III -22). For small ^ — - < 0.2} ; hence, 

Eq. (Ill -21) may be approximated by 















CL 



■ 



Ji 



(iv-1) 



With Eq. (IV-1), one can predict the release curve to be para- 
bolic in shape when f is plotted vs time to the one -half power; the 
curvature depends upon b. The curve consists of a straight-line com- 
ponent and a parabolic component, which is the evaporation correction. 

Figure IV. 1 is the graphical presentation of Eq. (IV-1) for those 
values of f that do not differ by more than approximately 0.001 from 
the more accurate expression for f. 

By assuming that b (the rate of boundary movement) has been 
determined, a corrected value of the fractional release may be ob- 
tained by subtracting the quantity .from the actual fractional 

release. This corrected value of f may then be used with Eq.(III-22) 
to determine D. 



E. Slab Case --Moving Boundary 



1 1/2 

For large values of — (3r , the fractional release is approxi- 

mated by (3T--the percent of the body that has sublimed. 

q t l/2 /g r l/2 \ 

For small values of t— ^ < 0.2 J , a simplified expres- 
sion for the fractional release may be developed in a manner similar 
to that for the spherical case. This is written as 
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6 




T 



1/2 



(1) 


P = 0 (no evaporation) 


(4) 


P = 200 


U) 


p = 20 


(5) 


P = 500 


(3) 


p = 100 


(6) 


p = 1000 




MU .32931 



Fig. IV. 1. Simplified expression for fractional release 
for the sphere. 
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i<&) 





(IV-Z) 



Hence, D may be determined in a manner similar to that for the 
spherical case. For small periods of time, the correction factor 
^ [>'!(¥)] is negligible; thus a straight-line plot may be ex- 
pected; and over longer time intervals the plot may be expected to 
show an increasing slope. 

F. Some Numerical Considerations 

Using the values of Shaked for his sample # 108 at Z040°C, one 
obtains 

D = (o. (o * I0' ,Z cm*/ sec, } 
t' /z = 150 (sec/ /z (~t = 2.Z5 * 10* 

and (£- ~ 0.2 5 * 

Therefore, T~ ~2 58 x /0 < ° (Z ^ = A 5 V v 10 ' 



The expected fractional release, with the use of Eq. (III-ZZ), would be 

■P-£i r' 4 - = 5.2wo' s - 

Shaked estimated that the value of b he encountered during the 
anneal of this sample was 0.3 p/h. Because Shaked's sample was an- 
nealed in a Knud sen -type cell, this evaporation rate was considerably 
lower than would have occurred if the sample were completely exposed 
to vacuum. 

Thus, 

| 3 - %-= 3 / 6 , 



and, if the value of f calculated above is substituted into Eq. (IV- 1), 
one has 
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^ 7 r 




5.2 X 10 ] 



with D representing a modified value of the diffusion constant 

D~ / x i0~ ' Z c ' 

which is 6 1% less than that calculated by neglecting evaporation. 

If the sample had been completely exposed to vacuum at this 
temperature, the evaporation rate, b, estimated from thermodynamic 
data, would be about 15 0 ji/h, from which 

jS = I5Q,000 ■ 

4 

For an anneal time of 2.Z5 X 10 sec and an actual D of 
6.6 X 10 ^ cm^/sec, it follows that 




Using these data in Eq. (III-2 1), one predicts a fractional release of 
f= 0.757. Now, if the experiment had been performed in this way, 
and if the effects of evaporation were ignored, these values of f, t, 
and a used in the nonevaporation Eq. (III-22) would yield a value for 
D of 

V- 13.9 X I0~ 3 CM 3¥ec / . 

which is approximately 21,000 times as large as the actual diffusion 
coefficient for this temperature. 

These calculations illustrate the importance of either minimiz- 
ing evaporation experimentally, or of taking it into account analyti- 
cally whenever one obtains diffusion data at high temperatures. 
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V. CONCLUSIONS 

An equation has been developed that relates the fractional re- 
lease, f, of a tracer element in a geometrical body under conditions 
of evaporation of the body and diffusion of the tracer through the mov- 
ing boundary of the body, with these conditions occurring simultane- 
ously. This may be better expressed as the fractional release of the 
tracer element through a moving boundary. 

For small values of the quantity |3, the release curve is approxi- 
mately a straight line whose slope is slightly greater than that for the 
stationary -boundary case. As (3 becomes larger in magnitude, the 
slope of the curve increases with increase in time and approaches a 
limiting value (see Fig. III. 2) for the values of (3 and t considered. 

It can be seen that to completely ignore the evaporation of the 

body may lead to an incorrect value for the diffusion constant. The 

1/2 

curve of f vs (t) ' can appear to be a straight line when the bound- 
ary motion is not zero, but the slope of this line is not the same as 
that for the ze ro -boundary -motion line; hence, different values of D 
would be obtained unless a correction were applied to the moving - 
boundary case. 

Simplified expressions have been developed that, for small 

l/2 

values of the variable (3r /2, may be used to make a more accurate 

determination of the diffusion constant, D. With these expressions 
one can predict a release curve whose slope, when f is plotted 
against the one -half power of time, will increase with time. The cor- 
rection made to the fractional -release equation for the stationary - 
boundary case is independent of D, and by subtracting this correction 
from the actual release, one should obtain a line of constant slope. 

It also should be possible to fit experimentally obtained 
fractional - release curves to the curves obtained in Sections II and III 
of this report, and thus to determine both the rate of boundary motion, 
b, and the diffusion constant, D. 

Experimental verification of the simplified and the approximate 
relations for the fractional release remains to be demonstrated. 
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APPENDIX 

A. "Conventional 11 Analysis of the Slab Problem 

Solutions to the heat-conduction equation for a slab -a<x<a, 
with a constant initial temperature Tq, and with T = 0 at x = ±a are 
given as^ 



-nvh 



and 




+eHc Kt^r j 



(A-2) 



The average temperature in the slab is given as 




■£, ex p t 



/C(3.n-f'l) 2 7T zr t 

^iOL- 



(A-3) 



and 




> 1-1 




(A -4) 



Using an analogy to the heat-conduction solution, it will be shown that 
the concentration of a tracer element C(x, t) for any time t and any 
position x, where -a<x<a, the initial condition C(x, 0) = C^, where 
Cq is a constant, is 

C (*t) = W*£$S- **pt ^P^ ]cOs(^^ (A-5) 

and 



»>■ c. - . ».« 
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Here, D is the diffusion coefficient for the particular element being 
studied. The corresponding fractional release, f, defined as the ratio 
of the decrease in average tracer concentration to the initial average 
concentration, is 







and 




( A-7) 



( A-8) 



A. 1. Tracer Concentration — Exponential Form 



By using the method of separation of variables, it may be 

g 

shown that a solution for the diffusion equation 





c) C (x it) — ~p (xfi) 


( A-9) 


with 


C(*,o)= Co 


(A-9a) 


and 


C - 0 J 


( A-9b) 


may be 


expressed as C (x, t) = F (x) G (t) , 
n ' n ' n 


(A-10) 


whe re 


F (x) = A cos k x + B sin k x 
n n n n n 


(A-lOa) 


and 


exp(-k»T>t) • 


(A-lOb) 



Applying the boundary conditions (A-9b), yields B^= 0 and 

k = r^(2nH), where n = 0, 1, 2, • • • . Applying the initial condition 
n c, a 

(A-9a) yields a = C^, 



so Eq. (A-10) becomes 
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C„(x,t)= C 0 A n cosk n x exp(-L z Vt) 



(A -ll) 



The coefficient may be evaluated by setting t = 0 and using the 



where f(x) = 1. (A-13) 



property of a Fourier series 7 that 

a n --j-j-f(x)cos-^-d x i 

a. 

thus> A n - dx * 

-a 

This may be integrated directly to yield 

a = » n = 0, 1, 2 • • ■ . 

n 7 T[2n+i) 

Therefore, the substitution of Eq. (A-14) into Eq. (A-ll) gives 

e ^L ice- r 2a 

Any sum of solutions is also a solution, so the final result is 

m- ■ 



(* -IZ) 



(A-14) 



(A- 15) 



(A- 16) 



For long periods of time, the tracer concentration should become 

By noting that the exponent in each term approaches 



zero as t 



OO . 
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-oo as t -► oo, it is clear that the tracer concentration will indeed ap- 
proach zero as a limit. 

A. 2. Fractional Release — Exponential Form 

The fractional release, f, may be determined from 

X- QjOLQJtl (a- i7) 

l - 9o 

where Q(t) is the total number of tracer atoms per unit area contained 
in the solid at any time t, and Qq = 2aCQ is the total number of tracer 
atoms per unit area contained in the solid at time t = 0: 





s&LL'o 




( A- 18) 




For t -> oo, the summation term approaches 0 as a limit; hence, the 
fractional release approaches 1 as a limit. For t=0, Eq. (A-18) be- 
come s 



-f 



= l-8_ p / 

TT Z Oio fcn+0' 



(A- 19) 



is shown 7 to be equal to 



The sum + + £ 1 

l2 32 52 n _Q (2n+l)^ 

it 2 8 it 2 

-g-, so f = 1 — X = 1-1 = 0 at time t = 0, which is as expected. 

7T 
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Thus, it has been demonstrated how Eqs. (A-5) and (A-7) were 
obtained, and it has further been shown that the tracer concentration 
and the fractional release approach the proper limits for large periods 
of time, and at time t = 0. 

It will now be shown, with the use of La Place transforms, that 

Eqs. (A-6) and (A-8) are the solutions of Eqs. (A-9) and (A-17), re - 
7 

spec tively . 



A, 3 . Tracer Concentration — Error - Function Form 
Let 9C(x,t) — 'T'l c>y (x,~t) 



with d (X, 0) - C 0 


( A-9a) 


and C &/£) — 0 j 


( A-9b) 


by symmetry, j — /C * 


( A-20) 


Now the La Place transform of C(x, t) = C(x, p); 




C(x,p)= f e' pt C(xt)<fr ■ 

J o 


( A-2 1) 


Making the transformations of Eqs. (A-9), (A-9b), 
tained 


and (A-20), one ob- 


pC(*,p)-C^T?'£^ Pi, 


( A-22) 


with Q (±a / p) — 0 


(A-ZZa) 



= o - 



and 



~bC (0, p) 

^x 



( A-22b) 
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Defining q‘"=p/D and by substituting in Eq. (A-ZZ), one obtains 

ycjx. p) (A-24) 

Solving Eq. (A-24) for C(x, p), one has 

C^p) = Aef -t- T3e~f-i- yr ■ < A - 25 > 

Upon applying Eq. (A-22a) and Eq. (A-22b) to Eq. (A-25), it follows 
that 




Ca (ef^eV) 

P (eTi-e'T) 




(A-Z6) 



Multiplying the first term top and bottom by 



e -qa 



one sees that 




( A-Z7) 



and upon expanding 



( 1 + e 



-Zqa. - 1 



by the binomial theorem, one finds 



/ + e Z T ] = / - € z f+ (e' 2< f/V ■ ■ ■ 

( A-Z8) 

= i(-/)Te-* n r. 

yi'o 



By substituting (A-28) into Eq. (A-27), it follows that 
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Taking the inverse transforms, one obtains 




( A-30) 



For long periods of time, the tracer concentration should be- 
come zero as t oo # If t is allowed to approach infinity in Eq. (A-30), 
the erfc — 1,0, and Eq„ (A-30) becomes 



&e(x,t)=c o -c 0 z(-ir(2) 

= C 0 -C 0 (2-2+2-2 + -)’ 



( A-3 1) 



which appears to be nonconve rgent and equals or -C^, depending 

upon whether the last term of the expansion is negative or positive. 
This is a periodic function, however; and if the function f(x) is defined 
as 



and 



f(x) = 2 


if 


0 « x < tr , 


(A-32a) 


f(x) = 0 


if 


-TT < x < 0 , 


( A-32b) 



then the coefficients of the Fourier series for f(x) are 

& 0 ~4ffo d* +-frf2dx -2 ’ < A - 33 *> 

-7T 

dx nx] o = 0 > (A- 33 b) 

and b n ~d-j o 2sir) nx dx = : ^—Cos nxjj= ^.(l-cosnir) > (a-33c) 

and the Fourier series for f(x) becomes 



oo \ 

-f(x) = ~^a 0 + 2 (&» Cos nx -hb„sm nx ) 

\ . d f sm x ,siy) 3x j. s m 5x 

= I + 1F(— r~ + ~3 5 
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9 

On page 180 of Sokolonikoff and Redheffer, it may be seen that f(x) 
is well-defined on the closed interval -tt ^ x ^ tt except for the simple 
discontinuity at x = 0, has a finite number of maxima and minima, 
and is defined for other values of x by the periodicity condition 
f (x -1- 2 tt ) = f (x) ; therefore, the Dirichlet conditions hold for f(x), and 
the Fourier series for f(x) converges to 

= +-f(x')] (A - 35) 

for every value of x. 

In particular, at x = 0, 

-C(o) = z i ~f(0 + -f 0 0')] = z( Z+0 ) = 1 ' <A ‘ 36) 

Because the coefficient of the second term in Eq. (A-31) is shown to 
converge to the value 1, the tracer concentration will, in fact, ap- 
proach zero as a limit as t -► <*> . 

A. 4. Fractional Release — Error -Function Form 

Calculating the fractional release from Eq. (A- 17), one obtains 

r 2aC 0 ~ / C(x,-i)dx 



a 
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which, with the use of Eq. (C-ll), becomes 

M' 



■f- 



( A- 38) 



- 2($)' i {z<-,r[ u ^ 

=4 $f[^c o ] . 

By letting n + i = m, the last two terms may be combined, and by 
noting, from Eq. (C-7), that ierfc 0 = tt , one finds Eq. (A-38) to be 






( A-39) 



. 1/2 



For t -»■ 0, the terms in Eq. (A-39) involving ierfc k/t 1 ' become zero 

thus, only the first term — ^ 7 — remains. For t = 0, this 

it 1/2 V a 2 / 

term and also the fractional release are zero. 



For large t, as t -► oo, the terms involving the 
na/( Dt) ierfc 0 

tions Eq. (A-39) becomes 



ierfc na/( Dt) ierfc 0, where ierfc 0 = tt with these substitu- 






( A-40) 



By using the same analysis as that for the tracer concentration, it can 
be shown that the summation term converges to - 1 ; thus Eq. (A-40) is 
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which is indeterminate; therefore, it is necessary to refer to Eq. (A-37) 
to resolve the indeterminateness. /“/enC6^ 






+ 2jJr l T(0~]dy, 








( A-4 1) 



This previously was shown to be 1, so 



lir* ~P — / 



( A-42) 



It has now been shown how Eqs. (A-6) and (A-8) were obtained, and it 
has further been demonstrated that the tracer concentration and the 
fractional release approach the proper limits for large periods of 
time, and at time t=0. 



A. 5 . Correlation Between the Two Forms for Fractional Release 

Using the method of Booth, ^ one can show that Eqs. (A-7) and 
(A-8) are identical for values of 7T — ^ < 1. Here, 



r fe^/)V 2 Dt 



p_ i / n-c f- ( 2 *r ,j 

T - l-Trz£l5^ eX PL ^ 0- 



(A-7) 



for 



. 2 Dt 



< 1. 



Equation (A-8) becomes 

f - z 2 

.1/2 ^ a 2/ 



( A-8a) 
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Rewriting Eq. (A-7), one has 

h^Tr^Dt 



■A 1-mi.e 



Ho?- 



, _ (Zmffr^Dt 

ZuLve^r- 



YY)-t 



(A-43) 



The sum g e ~^ x - / + 2 £ £ 
rr\~-c>o rtf-i 



-m z X 



oo O 

/ "XU“ 

e du. Thus, 



( A-44) 



-xu?- 



— oo 

where x for this case is 




\+zie-^ = (%f ■ 



( A-45) 



(A-46) 



Now, if both sides of Eq. (A-46) are integrated over the limits 0 to x, 
the result is 



x+22^0- 



m-i 




which may be arranged 





(A-47) 



( A-48) 



Upon substituting Eq. (A-48) into Eq. (A-43), one finds, 
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= 1 * 4 . fiSPf- 

_ _S f_L [ firr^pt y^ Rn±Dt_ y l 1 

ir^lH i( Ha.' 2 ' / aa* ^,w\| 

_ I + Ji ( T>t) 'A £ <? JL. . 

00 4 n 2 

c—* d. / TT 

The 2_, — is shown as equal to — ; hence, 



( A-49) 



n= 1 n^ 



f- 



tt'A V a?- J 



'/z 



(A-50) 



For values of 



TT 2 Dt 



< 1, Eq. (A-8) reduces to Eq. (A-50); thus, for 



these values of 



TT 2 Dt 



, Eqs. (A-7) and (A-8) are equivalent. 



Solutions of the form of Eq. (A-7) converge slowly for small 

values of Dt/a , but solutions of the form (A-8) are rapidly conver- 

7 

gent for the same small values. This statement also applies to the 
equations for the tracer concentration, i. e. , to Eqs. ( A-5) and ( A-6). 
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B. 'Conventional M Analysis of the Sphere Problem 

Solutions to the heat -conduction equation for a sphere, 0 < r < a, 
with a constant initial temperature Tq, with T = 0 at r = a, and with T 
finite at r = 0, are given as^ 




(B-l) 



( B - la) 



(B-2) 



( B -2a) 



[See note following Eq. (B-6a).] 

The average temperature in the sphere is given as 




and 





(B-4) 



Using an analogy to the heat -conduction solution, one can show that 
the concentration of a tracer element C(r, t) for any time t and any 
position r, where 0 < r < a and the initial condition C(r, 0) = Cq, 
with C Q a constant and with C(0, t) finite, is 
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ant) = s/ 






e (0, t) = C 0 i(-O n exp (-^ t ) » 



(B-5) 



( B -5a) 



and 



C(^= C. 




.J5 LfP&h 

ni 0 c nc Z(C 



l)a- 

z (Dt? /z 



r 




(B-6) 




( B -6a) 



[it will be noted that Eq. (B-6a) differs from ( B -Za) - -other than 
in the change of notation- -by a factor of two in the summation term. 

It will subsequently be shown that Eq. (B-6a) is the proper expression 
for the tracer concentration at the center of the sphere.] 

The corresponding equations for the fractional release of the 
tracer element are 




(B-7) 



(B-8) 



B. 1. Tracer Concentration — Exponential Form 



The radial -diffusion equation 





,2_ 3C(r,t) 

r 



(B-9) 
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may be reduced to the linera -diffusion equation (A-9) by making the 
substitution u(r, t)=rC(r, t) 



bu_(r,i.) _ _ ( r+) 

” 3 t -P ’ 






( B - 10) 



with the following boundary and initial conditions 

u(r, 0) = rC Q , 



and 



u (a, t) = 0 
u (0, t ) = 0 



(B-lOa) 

(B-lOb) 

(B-lOc) 



A solution for Eq. (B-10) may be obtained by using the method of 

Q 

separation of variables 0 



where 



and 



U (r, t) = F (r) G (t) , 
n ' n n 

F (r) = A cos k r + B sin k r , 
n n n n n 



G (t) = a e 



-k 2 Dt 
n 



n 



( B - 1 1) 
(B-lla) 

(B-llb) 



Applying the boundary condition (B-lOc) yields A =0. The boundary 



condition ( B - 10b) gives k =^-, whereas Eq. ( B - 10a) gives a=C f ,, 

n a u 

(n = 1, 2, 3 • • * ), and Eq. (B-ll) becomes 



n 



(J„(r,t) = C 0 B„ sin k„r exp (-kfvt) 

= C a B n smfr exp f-^j. 



(B- 12) 



The coefficient B^ may be evaluated by setting t = 0 and by using the 



property of a Fourier series that" 



V) J 



JT 



(B-13) 
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Thus, 




sin 



MUX <i r , 
a u 



where f(r) = r. (B-14) 



This may be integrated by parts to give 




irn 




( n = 1, 2, 3, • • • )• 



( B - 15) 



The substitution of Eq. (B-15) into Eq. (B-12) results in 

U n ( r /'t) ~ ~~ 2O-Co(~0 eXp S/n ~~ • (B- 16) 

Any sum of solutions is also a solution; so the final result, upon using 
C(r, t) = J u ( r, t), is 

am) ( 4 V [-*%&]*'■" 

and 

where the relation 



(B-17) 



( B - 17a) 



lim sin x 
x — *0 x 



has been used. 

For long periods of time, the tracer concentration should be- 
come zero as t -*• oo. By noting that the exponent in each term in both 
Eqs. (B-17) and (B-17a) approaches -oo as t °o, it is clear that 
the tracer concentration will indeed approach zero as a limit. 
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B.2. 



Fractional Release — Exponential Form 

The fractional release may be determined from 



-f= 



_ GL-QW 

do 



(A- 17) 



where Q(t) is the total number of tracer atoms contained in the solid 

4 3 

at any time t, and Cq is the total number of tracer atoms 

contained in the solid at time t = 0. Thus, 



f= 



% rra 3 C 0 - J o ^-n-r z C(n i: )^ r 



_ r/ 3 



I- 



V / 3 TT- a 3 Co 

3 (. ZaCte C-)? 






TT 



^ V ~ Y 



rfVty. 



r son dr (B _ 18) 



_ 6. _L 



7T~?j='l 



J? r!* ex P (- 



n z rr*Dt ' 

03 J 



For t ->oo, the summation term approaches zero as a limit; hence, 
the fractional release approaches 1 as a limit. For t = 0, Eq. (B-18) 
become s 






Ill 

The sum + + + 

l 2 2 2 3 2 



(B-19) 



00 i g ^ 

= V is shown 7 to be equal to ; 

*—> "> A 



n= 1 m 



so, 




1 = 0 at time t = 0, which is as expected. 



TT 

Thus, it has been demonstrated how Eq. (B-5) and (B-7) were 
obtained. It has further been shown that the tracer concentration and 
the fractional release approach the proper limits for long periods of 
time and at time t = 0. 

7 

It will now be shown, by using Ea Place transforms, that 
Eqs. (B-6) and (B-8) are the solutions of Eqs. (B-9) and (A-17), re- 
spectively. 
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B.3. Tracer Concentration — Error -Function Form 




We have 






XXt) _ j, px(nH + 2 
~5t' L) i dr* r 


'S ^ 
pN 

i / 

v» 


(B-9) 


with C ( a, t ) = 0 , 




( B -9a) 


C ( r, 0) = C Q , 




( B -9b) 


and 0 9 by symmetry . 


( B -20) 


The La Place transform of C(r, t)=C(r, p) 


is 




C(r,p)^ Je~ pt C(r l t)dt ■ 

o 




( B-2 1) 



Making the transformations of Eqs. (B-9), (B-9a), and(B-20), one ob- 
tains 



pC(r,p)-C o = 




( B -22) 



where C(a, p) = 0 , 

and 2^ JCQJ?). — Q . 

dr 

Defining 

q 2 = P/D 



and substituting Eq. (B-23) in Eq. (B-22), one has 




Solving Eq. (B-24) for C(r, p), one obtains 

C(V,p) = -£ eV-f-Y-e'V+js- ■ 



( B -22a) 
( B -22b) 

(B-23) 

(B-24) 



(B-25) 



-53- 



Upon applying Eqs. (B-22a) and (B-22b) to Eq. (B-25), it follows that 




CqQ [e ^-er] Co . 

rp [e-f- eVJ P 



(B-26) 



Then, multiplying the first term of Eq. (B-26) top and bottom by 
e ^ a , one finds 



C(V, P )=- e^][l- e z f] 





( B-27) 



and expanding [l - e 2qa ] * by the binomial theorem, one gets 



]-e = I + e ^ a '+^.~ z T) z +(e z V : f-t- ' 

= i? e- 2n ? a • 

r i-o 



Substituting Eq. (B-28) into Eq. (B-27), one sees that 



( B-28) 




rp 








( B -29) 



Upon taking the inverse transform, one finds 



C(rt) = r 0 ^ f^r r (2»+lh-r r (Zn+px+r 

L\ r >L) C 0 r fe, (_ S+C 2(Dt/» erfc 2(K)* 



( B -30) 



and 

C(0,t)= r? 0 C(r,t) = Co-# > i B - 50a > 

10 

which is indeterminate. Applying L, 'Hospital 1 s rule to the second 
term to eliminate the indeterminateness, one sees that 



lim ( x ) - hyw f C x ) _ 

^(x) X-*a. g ' fx) 



• 44 
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Now, 

C(o,t)-~ Co- !i* | o [er/c%|^- e^%^fe r 




As an additional check on Eq. (B-30b), it can be shown that the 
tracer concentration approaches zero as a limit as t -►oo. 

The sum 



£ e~ hZx = 



l+2£e~ li * x 

n=i 



is, for small x, essentially J e XU du, 

-Too 

where J 6 = (gf . 

— QO 



(A-44) 



( A-45) 



Equation (B-30b) may be written as 




(B-31) 



Combining Eqs. (A-44) and (A-45), one finds that 






/ 



5 ^ _ _l(jt __ 

~ zUj ' a 



( B-32) 
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and by substituting Eq. (B-32) into Eq. (B-31), one obtains 




= Co - Co = 0 • 



B.4. Fractional -Release — Error-Function Form 



Calculating the fractional release from Eq. (A- 17), one has 




Now, Eq. (B-34) may be integrated by parts, by using Eq. (C-ll) and 
by letting u = r and 

dv= erfc (2 z(l±y£dr - 



which yields 




-'5 6- 



&(%*) [M Ler ^ c 1 e ] 



(B-35) 






By letting n + 1 = m, the ierfc terms may be combined; the 

Z Z 

i erfc terms will cancel except for an i erfc 0 term. Noting that 

the ierfc 0 = tt and i^erfc 0 = l/4, one finds that Eq. (B-35) be- 

come s 



T~ 6(§f^" /z +zSier^^~l - 5 ^) • 



( B-36) 



. 1/2 



For t -► 0, the terms in Eq. (B-36) involving ierfc k/t A/ “ become 
zero. For t=0, the remaining terms are zero; thus, the fractional 
release is zero. 



If the substitution of Eq. (C-7) is made for ierfc 
Eq. (B-36), the resulting expression becomes 



na 



VDt 



m 



f=6(f) /2 



T T' k +Zn ' k z exp (- -pr ) 



(B-37) 



For large t, as t oo , -*• 0; thus, Eq. (B-32) may be substituted 

into Eq. (B-37). Hence, 



linj X= 

iz-?o* f->oo 



[s (■§■/ (" r ' ■ zi (&§& e ^ 

,/z [^(rrDtf- l] j - 3 @)j 



+ TT 



— I/m 



[-/2|(n)erfc^+6(^)-3(i)] 



( B -38) 



— hr w 



5 ^) _l2 n?, (n)er ^ C 



nd. 

(D t)'/* 
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This is equivalent to oo -> oo, which is indeterminate. Referring to 
Eq. (B-34), one finds for large values of t that the difference be- 
tween the two complementary error functions may be approximated by 



Z Ifenilktr 



fen+l)cL-ir 



exp 



(Zn+l)' 



- 2]r_ r 6nt/)V 
' (TTDt)^ eX P L VDt 



which, upon substitution in Eq. (B-34), gives 

fcj = fife- Zex P 



— iim 



6 



t-**° 0^(TTl>t) ,/7 - J 0 (n= 



Applying Eq„ (B-3Z), one finds 



a 






= r 2 [2( nVt )‘ k - (”&>*] dr 






— — \ r z dr ~ — / • 

(3 3 V 34* 1 



Thus, l im £=1. ( B -39) 

t~>00 

It has now been shown how Eqs. (B-6) and (B-8) were obtained, and it 
has further been demonstrated that the tracer concentration and the 
fractional release approach the proper limits for large periods of 
time and at time t = 0. 

2 ^ 

Booth has shown that for values of — — - — < 1, Eqs. (B-7) and 

6 aZ 

(B-8) are identical. 



-58- 



Solutions of the form (B-7) converge slowly for small values of 

Dt/a , but solutions of the form (B-8) are rapidly convergent for the 

7 

same small values. This statement also applies to the equations for 
the tracer concentrations, i. e. , to Eqs. (B-5) and (B-6). 
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C. Properties of the Error Function and Related Functions 

7 

The following definitions are made: 

x 2 



erf x = f e Z dz , 

*frr o 

erf oo = 1 , 
erf ( -x) = -erf x , 



(C-l) 



(C-2) 



( C -3) 



2 r 3 -z 2 

erfc x = 1 - erf x = j e dz , 

\TtT x 



(C-4) 



ierfc x = erfc z dz , 
x 



(C-5) 



.n r 

i erfc x 



00 

= / I "" 1 



erfc z dz , 



(C -6) 



where 



.0 , - 

i erfc x = erfc x , 



ierfc x = — — e X - x X erfc x , 



nTtt 



( C -6a) 
(C-7) 



and 



A 

i^erfc x = — [erfc x - Zx ierfc x] 
4 



The general recursion formula is 



2n i n erfc x = i n-2 erfc x - 2x i n *erfc x 
If both sides of Eq. (C-6) are differentiated, 

i n erfc y dy 
z 

dz 



.n- 1 



d .n+1 c d 

l erfc z = — — 

dx dx 



.n , 

= -l erfc z 



dx 



(C-8) 



(C-9) 



(C - 10) 
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Integrating both sides of (C-10) and rearranging yields 



I 



• n C 

i erfc z 



dx = - 



1 

dz 

dx 



. n+ 1 

l 



erfc 



z 



(C-ll) 



provided — is a constant and the constant of integration is neglected; 
dx 

then (C-ll) is the indefinite integral of i erfc z dx, where z is a 
function of x. 

The values of i n erfc(-x) will now be studied to obtain an inver- 
sion formula, so that the tabulated values of the positive argument 
may be used. 



Noting 


erf ( -x) = -erf x 




(C-3) 


and 


erfc x = 1 - erf x , 




(C-4) 


one obtains 


erfc ( -x) = 1 - erf ( -x) = 1 + erf x 






= 1 + ( 1 - erfc x) 




(C - 12) 




= 2 - erfc x . 






Thus, 


f 1 - x ^ f 

lerfc x = e - x erfc x , 

n/~tT 


(C-7) 


and 


1 -( -x) ^ 

ierfc ( -x) = e - ( -x)e rfc ( -x) 

\TtT 






1 -X 2 

= e + xerfc ( -x) 

n/~tt 

1 -X 2 

= e + x ( 2 - e rfc x) 

\Ttt 




(C - 13) 




, , r 1 _x2 f 

J 








= 2x + ierfc x . 






So, 


i 2 erfc x = — [erfc x - 2x ierfc 
4 


x] , 


(C-8) 
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and 



inerfc ( -x) = — [erfc ( -x) -2 ( -x) ierfc ( -x)] 
4 

= — [ ( 2 - erfcx) + 2x ( 2x + ie rfc x) ] 

4 

= — - — (erfc x - 2x ierfc x) + x^ 

2 4 

= — + x^ - inerfc x . 

2 



(C - 14) 



By using the general recursion formula (C-9) and the relations just 
developed, it may be shown that succeeding terms are 



.3 , , . 1 , 1 3 , .3 , 

i erfc ( -x) = — x+— x +i erfc x , 

2 3 



(C - 15) 



.4 , , , 1.1^2,14.4. 

i erfc ( -x) = — + — x + — x - i erfc x , 

16 4 12 



(C - 16) 



and 



.5 , , x 1 , 1 3 ^ 1 5 ^ .5 , 

i erfc ( -x) = — - x + — x + — x + 1 erfc x 

' 16 12 60 



( C - 17) 



From the foregoing, a general inversion formula for i n erfc(-x) may 
be obtained: 



i n erfc ( -x) = ( -l) n+ ^ i n erfc x + + — 

n! 2(n-2) ! 

n-4 

, x 



(C - 18) 



2x 



n-m 



16( n-4) ! 



(4*8 '12* •• 2m) (n-m) ! 



whe re 


n - m = 1 


for 


n = 3, 


5, 


7, • • •, 


and 


n - m = 0 


for 


n = 2, 


4, 


6, • • • 
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D. Solution of the Slab Problem for a Moving Boundary 

DM- Concentration of the Tracer 

An expression for the fractional release of a tracer element 

„ 5 

through the moving boundary of a slab will be obtained. Chambre 
has derived a general solution for the tracer concentration for a slab 
-(a - bt) < x < a - bt and 0 < t < — , where b is the rate of movement 



of the boundary ^ 

C(x,t)-2^ Wz fy)eyp[-Ml d J (D1) 

+ 2(-l)"JfCs) exp [•# (s +na )J S M 



s (x,±;j) = exp [-(£$] + 




(D-la) 


This equation is the solution of the diffusion equation 




ciCCxyt) _ -p, ^c(xX) 
bt " ^ 




(D-2) 


subject to the restricting conditions 






C6m)=T 6<) for 


0 


( D-2a) 


C)C(Ot) _ /\ for 

~Jx U 


0< t<£. 


( D-2b) 


and 0 for 


0 <t <•%•• 


( D -2c) 



For this problem, the initial condition (D-2a) shall be assumed to be 
a constant, C^; thus, 

C(k o) = Co Ior o<lxl<a- (D-3) 



If this substitution and S(x, t; £ + 2na) are made into Eq. (D-l), one 
has 
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a. 



+ M-/) n expj 



m-\ 



( e *P (- + exp [~ ] • ( D - 4 1 



This may be simplified somewhat by defining 

y = £ + Zna 



and 



upon which Eq. (D-4) becomes 



x 



- £ 



2 \TDt 



(D-5) 

(D-6) 



(D-7) 



X—CL* 

, A?5t)K 

c (x V = rffem (-^y e *p (- 

2^)'^ 

The first term may now be integrated directly, and if exponentials are 
combined, one has 

cyy = % f erf [$$*] - ^[m 4 . 



-f 



(Z/i+f) 

X. I 

af2»- 1) 



1 } 



,d - 8 > 



+ 






The exponentials may be simplified, 
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= ^jjj -( x ~ t 2r l yt)'J _ (x+hbt)~^ 

~ vit k -( xtZnht ^ ~ ^ 



( D-9a) 





t 



} 



_ / 

yz>t 

~w 



jjpt-2^(x-2hbt) + (x~2nkt) -hx z ~(x-2hbt) z 'j 
• [(j+(x-2hbt)] + %bt(x-nht)^j 



( D-9b) 







Simple substitutions similar to Eq. (D-6) allow the last two terms of 
Eq. (D-10) to be integrated: 



-65 - 



c(x,t)= [er-f^r x - erf^n] 



— €r~ f a(2n-l)-(x^2nht) 
2 [m^ 



(D- 1 1) 



_ r <z(2n-i)+(x-2nbt) l 

2 CDt) /x J 

Using Eq. (C-3) and rearranging terms, one has 

c (*, t) = -jf [e ^^ 5 + erf w>] 

+ 2^S( L i)" e X p[-^ (m-M-x)]je rf 

-e rf ? ^M>y £r* j (d-ui 

With the use of Eq. (C-4), i„ e. , erfc x = 1 - erf x, one finds 

eft 9 = C, - [ er * J#- +er ^ &t] 

+ ^§(-O n expl:& (m - ntt -{l[ er ' f ' c M z($f£~ Y 

- erfc 2w £ a ;te *- 



-f 



/C ^7-/ ^ 

-erf c 2 &d£p$ 4 *~l. 
enc Z(Dt)'& J 

For t = 0, because erfc oo = 0, Eq. (D-13) becomes 



(D-13) 



C(x, 0) = C 



0 
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For t = a/b, which corresponds to total evaporation, the tracer con- 
centration should become zero. Hence, 



C(*, a 



F 



)=c 



-A. 

2 







a~-x 



HW* + er ^ c 2 (- 




(D - 14) 




Since | x | < a - bt at t = a/b, x in Eq. ( D - 14) must be zero. Thus, 




with the use of Eq. (C-12), 




( D- 15) 



OO 

It was previously shown in Eq, (A-4 1) that Z ^ (-1) = 1; similarly, 

oo n=0 

it can be shown that Z ^ (-l) n =-l. With this substitution, 

n= 1 



Eq. ( D- 15) is 



CK )- = 0 • 



(D-16) 
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If, in Eq ( D- 13), b = 0, which co rresponds to no evaporation, then 

c (x£) =C,-^[er-fc -f- erfc %=& ] 






Co £(-lf[er-fc & - er r r (z^i)a-y 
Z ft, 7 L c,rc z (W& ertc z(d t)& 



( D- 17) 



+ e M 2n ;¥t x - erfc C 2ri+ P a+X 1 . 

+ TC 2(6t)^ ^ nc Z{ W)'/z J 



If 2n - 1-Zm+l, Eq„(D-17) becomes 




which may be recognized as Eq. (A-6). 



D,2. Fractional Release of the Tracer 

The fractional release may be computed from Eq. (A-17), as 
follows : 




<2o 



Xc i-kb) 

3a C 0 ~j C(x,t)dx 



2clC ( 
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using the symmetry relation, one obtains 



(&-&) 



(D-19) 



With the substitution of Eq (D-13) for C (x, t), it follows that 
.(a-b€) J /fci-fct) 



t D-20) 



l-aF 0 jS° dx 

- ^ f^l(-iTexp[-^-^[erfc 

_ p Ca-bt)+a-hX \ (iv 

Cr+C — J«<- 

With the use of Eq, (C-ll) in the second integral, Eq, (D-20) 
become s 

f ^ +^ fe f' terfC +t6 ^ W +t£rfc ^-' MC W‘ 

-^(-fexp H h m-yt)]fexp(-^rfc *t$x — 
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The remaining integrals in the expression (D-21) for f may be 
evaluated by integrating by parts, as follows. Let 




The result may be written as 




£ ( Dty* 



(D-22) 



To integrate the second integral by parts, let 




and 



dv- e rfe 2n ^ a ~^ ± — - dx 

av ~ erTC 



then 




and 
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The result is 




( D-23) 



Upon applying the limits to Eqs. (D-22) and (D-23) and inserting the re- 
sults into Eq. (D-21), one obtains 




After eliminating equal and opposite terms and combining exponen- 
tials, it follows that 
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( D-25) 




At t = 0, each term of Eq. (D-25) becomes zero; thus, the fractional 
release is zero. At time, t = a/b, which corresponds to complete 
evaporation of the solid, the fractional release should be 1, as shown: 



f = / + 2 [‘ er/c 



z(^) l/z 




■L^er-fc 



a 



z(W 



+ L^er-fc 



CL 



zwr* 



- L”e 



ir-fo. 



— a. 



zm 



& 



(D-26) 




If, in Eq. (D-25), b is set equal to zero, which corresponds to no 
evaporation, then 






_ (m 

(Z- 



J/z 



Ler-fc 0- Le r/c^gs. ] + (-§") ^ 



(D-27) 



tier-fcH^ -ierk^ 
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by letting n - 1 = p and n + 1 = m in Eq ( D-27), one obtains 

~f - (j?f Z \} er ^ c 0 ~ ier ^ c 7W* + cer ^ c 

i-^2('0 P+Z Cerfc-^z J 

- (^}' /z j2i er-fc 0 + 2%(->) oer-f C j^ /z 

■f 1 ier ^ c joDp- t ff~ ^ ie r ^ c jc2) 11 *' J 

-1/2 

Noting that ierfc 0 = n 7 , and combining the m and p terms, one 



(D-28) 



finds 



f = 2(w) /Z [^ 2£(-lf Cer-fc 

which may be recognized as Eq, (A-8). 
D„3. Conversion to Dimensionless Form 



D CL 



( D-29) 



If one defines 



and 



so that 



_ ba Tfc _ bt 



D a > 



T = 


Dt 

a^ 


( D-30) 


P = 


ba 
~D ’ 


( D- 3 1) 


bt 

a. 


) 


0 < 1 ^ ( D -32) 


t = o, 


and 


(3t = 1 corresponds to com- 



plete evaporation, then Eq, ( D - 1 3 ) becomes 

C(i, t) = C 0 - (erfz 2$ + erfc 2^) -/ 
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/ - / 






c~o r 

I C ^ <~*V » * , . . ^ V~\ O • 

r-~-^jrij exp / "/-< 

/C n~i ‘ ' 1 



T 



nKwHt)} [e rfc (D .„, 

- €v/c ^O^ldzi 1 

W~ r U 7z J 

IV' * /r\ , A | V j P -<A ( l-/3 , c)- l -f-aT 

c( - '-''a-h i: Xr 7 ^ 

- £>X 2r L Q-(3 r L)+\+% : 1 

2 Z J ' 

Equation (D-25) becomes 

P^pZ + Z^Lerfc-gfa-Cerfc %^) 

-I 



+ z 



, OO y^_/ 

^ ,/z jZ(-l) (Znpz U J ex p [- njtfn-iXl ~{ft) j[l m e 'rfc ^dlL^ 




(D-34) 



+ Z /z <2(->) (2r>/SC >h ) expf~H/6(n-hlXl-^) l]i rn er(c 1 

r, ,/>■■-. , I - I I fC L, "' 

- L^e rfc ^iZZtr^ZhL 1 , 

< c '" J 



and Eq. (D-29) becomes 



f- 



_ q r r'/z T~-'h 

~ xC - j/i 



^ ^(- if Cerfc 



jz. i 
r 7 - ) 



( D -25 ) 



D.4. Approximate Solution for the Fractional Release 



If Eq. (D-34) is expanded to the first few terms in n, then 



-f : =pZ + Z' u (ier-fc-^7i -oerfcZZ% ) 



U v (~yC - y ( V-/ W / ; rr'f~ v CHX 2- r fc' 12 - 



r/l-j 
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-f 



Now, if the restrictions — < 0.01 and (3 r < 0.5 are made, then 

a^ 

Eq. (D-36) may be closely approximated by 




(D-37) 



Inversion formulae for i m erfc ( -x) have been developed in Appendix C. 

If these are applied to Eq. (D-37) and expanded in terms of m, it can 

/o t 1/2\x 

be seen, if sufficient terms are taken, that terms of ( t- — ) cancel 



and only the i m erfc 



(3t 1/2 



terms are left. Thereupon the approxi- 



mate expression for the fractional release, f, becomes 






(D- 38) 
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E. Solution of the Sphere Problem for a Moving Boundary 
E.l. Concentration of the Tracer 

In this section an expression for the fractional release of a 

tracer element through the moving boundary of a sphere is obtained, 
x 5 

Chambre has derived a general solution for the tracer concentration 
for a sphere 0 < r < a - bt and 0 < t < a/b, where b is the rate of 
movement of the boundary 

C (n -t) 

( E - 1) 
(E-la) 



(E-Z) 



( E -2a) 
( E -2b) 
( E -2c) 

For this problem the initial condition, (E-2a), shall be assumed to be 
a constant, C^; here, 

C(r, 0) = C Q 




This equation is the solution to the diffusion equation 

Xfr.t) -p r ^CO.t) , ,2 XXX) 7 

at f dr J 

subject to the restricting conditions 

C(r, 0) = f ( r ) for 0 < r < a , 

C ( 0, t) is finite , 

and C ( a - bt, t) = 0 for 0 < t < a/b . 



for 



0 < r < a . 



(E-3) 
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If this substitution and S(r, t; £ + Zna) are made into Eq. (E-l), then 




The first term may be integrated by parts, by letting 

u -~5 



if the second term is rewritten, then 




The exponentials may be simplified, as follows: 



[r- o^ha)] / 

~£T~HDt 







- ¥mr+ ¥/?*%*+ r z ] 



=7f^lf-^^tf-2 S (r-Zna + 2M) +f (E _ 6a) 

+ T z - (r- 2na.+ 2nbt) z - Hnar 

- xy^T'/f r-2ra.+2fibt)-j] If- tnrfa-Lt) + L fn z (o--btf 

-f— Hnar- r 2 ] 

= j(, r-Zna. +2rt U)-j] ^-^r- [~ ^nar-f-Vnlotr + % 2 o8 

- 82dt+ Hnbtfbf) +Hnar- 2 ot ] 

^ -J 

and - -L-lrt+ZrCj+Znayfrz^-Mtj] 

— 7j^8 Ef^ - V/?aj - 8hbt^2 l y 

— —L- (f 2 +2f( 'T- t ’2ncL-2nkt)+ Hn ar + J iG x d z -hG z ~] 

~ fDt ^^nbt) +( r +Z^-2r>bi) Z 

— (n-2hd~ 2nbt) z + Hnar + L lr?o}+'r z ] 

— T^jyjr [f+C r+2na - 2n bt)J Z ~ 7^ [r z -*-¥hr (a- lot) + 

- Hnar - L fn z a z -r z ~) 

~7f^\j + ( r +2r)G.-2nh±) ^nor-Hubtr +Hr2a2 

- 8 r2rx\ot+ ^M(nbt) - ^nar-^n 2 ^ 
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Integrating the second term of Eq. ( E-5) with the aid of Eq. (C-ll), 
and also substituting the results of Eqs, (E-6a) and ( E-6b), one finds 




r-a 

2(Dt)‘ k 



+ er-fc 



a+r 

Z(^_ 








The remaining integrals may be solved as before by integrating by 
parts. Let . , 

U J 

and lef Jv = exp [- dj 



>d dv= 7 ?-^ dj , successively; then 

C(r,t)= ^[erfc^+erfrj^yZ 

- [ierfc o er-fc ] 

+ er-fc r ~ZjZ 2 Z +a ~ ] 

(\a ^ 

2r 



-f 



Zexp[-#Mtrr)][&fc r f^ ~ ^ ^0^ ’ 



-79- 



-^gexp^Z-rJjferZc <lj 

-CL 



(£-5) 



_ c 



£^gexp[^(no.-ntt+^jer& C±2^^ 

Upon integration of the last two terms of (E-8), one has 

C(r,t)= 

- £*$ 0 *[ieirfc 



+ 



C_ 

Zr 



+£ 

T z 



(&£ expf-^m-nit-xE er fc r 

^£exp[-^-^r)][er-fc r+ -f/rtV^" ~ (£ ' 9) 

r r ^ ' r r +2na-2n)ot-^ l 

~<L 



, r+ZtorZnbT-oJ] 

-+enc~- ^ J 



- e r-fc 

n ~ l / & lJL 'c-ZYiaJ-2n^tboJ I 

cerft - Z(Dt)'/z J 



r+2ncL-2hbt-cC 1 

cerfc z(Dk yfc — J- 



By using Eq. (C-1Z), one obtains 



r r-Zno.+Znb't-CL , „ _r r-2no.+2nk>t~ho~ 

eric 2X5§fe — +erf c ETcZJ 17 *- 

_ o sirr.2no-2ftbt+a.-r , o r 2no.-2nPt-r-a. 

-Z-er+c +2- eric Z (Dtfe 



(E- 10a) 



er ^ C £<$& ' 



z(oty* 



and 



(E-lOb) 
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Using Eq„ (C-13), one has 

lerSr* r-2ha.+ 2jnjst-(X. _ ; r 

cc r rc ^ r-nbv/z terrc 



Z(Tt) ,/z 

of / 2 na.- 2 ,infct— )r + G~ 

-2[—2fUEJ^ 



_i_ ; ^ uSZ. 2r<<x-2nbt+0--r 
^ 6e ^c — z(Dt)y* 



¥* 



^ 2 fM 

= + terfc. 3 %$$}fc a " K - cer-fc *%$$£ * * 



(E - 11a) 



and 



terfc im*- 



(E-llb) 



The substitution of Eqs. (E-10) and (E-ll) give 

c(r,±)=(zm) 

+ (tfzp) H : ex P [- J $( r,a -- r,bt - r )] 

-gfj?exp[-&rn-»&-r)][er-fc 

+erf c 2™- ( ^-°-- r ] 



+ JJWp [~ (n&~nkt+r)Jle hfc 



-h€r-/2 






2. 
hbt~ 






y- 
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— expf-4 (no.-nj^t-tj i] [ - ie ^ Zna^z^kt+^-l: 

r n=i rL u JL zTopyfe 

— i€ffc 2ng-2n LT- a.- r 
2(0 ty/x 

gj exppgfa- 



( E - 12) 



4 .CJPtf * ~ 



r 



-not+r 



A1 r? p v T r 2 na -2h2t +g + r 

jjiccnc * 

— i &irX r 2na.- 2nJ?t~d+rl 

c 'zCotyf ^ — J* 

Upon cancelation and rearrangement, one has 

c (r,t) = c 0 - [erfc%&- erfc 



-t 



"zr"^, ex P £"$ {nornht- r)][erfc — §{(§ ) 2 r 

-i- an £% 2/7 (&-&£) ~CL-/~1 

^ rtC ~2fS t)'M J 

^^exp[-^(n a -nbt^epc^^^ 

, r ^r zn(a.-ht)-a.+r 1 

k /■ ei ^sm k J 

\ierfc0^. - Ler-fc 



C&Cl ( 



Dt 

0?! 






( E - 13) 



-C^^J^e X p[-<^(na-nli-r)][ierfc 



+ 



- ier-fc gn £0jvr= zr ] 
£^(^fSexp[-^-^r)l[cerfc Zh ^k^ r 

- cerfc *i 2 $^ } 



For t = 0, since the erfc oo = 0, Eq. (E-13) becomes C(r, 0) = C^. 
For t = a/b, which corresponds to total evaporation, the tracer con- 
centration becomes zero. (The proof will not be given here.) 

If, in Eq. (E-13), b = 0, which corresponds to no evaporation, 

then 

c(r,t)= c 0 - ir- [ er ^ w* " ~ 
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CnO. < F er r c (2h+l)a.-r , ^ fen-/)a-r 
2r h L erfC z(Dm ^ eHc zFDE)^ 



_ C 0 CL(J2tfi 
r \ as 



- - er{c 



(5 




^ — / p i~~Cr‘ ClS )T~ I 

2(oty^ 6 e r ^ c 



■ 



- ^(S/I ^ (2 ^r 

- if rJ~~ (zm+l)a.j-r + (zn-i)a+r 7 

r+c ZjD^ + Ler+C - 2(Dt)’^ f 

If 2n - 1 = 2m + 1, Eq. ( E - 14) become s 

CSt) = Co ~ 2 ^" \f r ^2(Dfyz ~ er ^2$ff* 



] 



- C„g. g r^r Szn+i)a.-r „ r _ (J w /) a+r - 

2r fcl? r+c 2{5m erhl Z(ix/z J 

— Co G-^ 1 j~p r fr - C^rV'l}_ a [7 r — , a^X. (km+l)d+S 

2s& 0 L^ rrc zmvr er ^ 1 ^f5m 



7 



_ CoO-fM.S 

~{d % J 






cesfc 



g±r 

2(Dt) 



'/a. 



*wi '■ 









/i - C 0 a. g [eSr ^ n Sa 
L ° 2r fc 0 L- c 2fDt)''* 

, [erfct^fj&f- ~ erfc- 



2. G m? 0 

'%)% F erfc %$^ - ier ^ 



C 0 a 



r 



-14) 



-83- 



+ 




_ i&r-j^fem+Oa+r 



(E-15) 



— 0 — Cjjx_ g r arf , (zn+Oa-r fen+Da+r-] 

r L en zWE)^ cn '- 2(0 t) ,/z J’ 



which may be recognized as Eq. (B-6). 

E.2. Fractional Release of the Tracer 

The fractional release may be computed from Eq. (A-17), 

O 0 -Q(t) 

f = — , (A-17) 



where O(t) is the total number of tracer atoms contained in the solid 

4 3 

at any time t, and Qg = — Tra Cg is the total number of tracer atoms 
contained in the solid at time t = 0. 



fa-it) 



Thus, ~r — 



n ¥/ 3 it^Co ~ l ^7rr z C(^ )dr _ 
~~ 1/3 Tra 3 Co 



(<x-yO 



l~m;y C(r ' t ) dr - (E - 16) 

■vi 
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~h 



ifc fr erfc - ce r/c $&]• dr 

“ \^0j k £ exp [-f -fcerfc 

- Ler-fc ? 2(b$A * C ] dr 
] r z (¥i ^f^exp[-f(na-M + r)][ie rfc 

-Cer-fc^^i^Jdr 



dX, 

__ 3 _ 

a?Cc 



(a-bt) 

r 

(a-bb) 



= i -§ l rZdr + ^ I r L erfc »- er ( c m^ d 



3 



O?j 0 
Ca-yt) 



2^ 
3 



M-bt) 



~h 



2(btf*- 

" ' -/• er/c 

'^exp^( na ~ n ^ +r )][er-fc 

^ + e ^ ^ r0?^ ] d 

w &t{ r (% &p rfc ~ |^~ r 

- terfe 



(E - 17) 






r 



r 



- # (wT0, ex P c - 7 ^tP 

- -— 

mf0 ^ §&-^ fc ^ Y r - 



+-2 

0: 



The first integral in Eq. ( E - 1 7 ) may be handled directly, whereas the 
second and last integrals may be handled by integrating by parts twice, 
as follows with the use of Eq (C- ll) 

IJL- r > dVf - ^ er ^ c z(^ 7z w ^ere n = 0 } l, 
du.-dir • and V, ~ 2(Dty z i nt erfc ’ 
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or 



dv^Cerfc^^dr, 



Hence 



. f = I - i? (0--^f^[Z( Dty>][(«.-I>i) i e rfejty* 

^a- bt)i.erfc^J\ 

~Z£- !?^ 4 J[" ^OtA’fer/c^ijj + 2fDt) & 6 ^erfc j 



+ L * erfc m* + (a -' kt) ^ rfc . 

, (a-bt) 

+ .Jr 2 exp[-#(n<i-» it- 1 -)] [erfil 

+ erfc^0^]dr 

- £ Jr 2 exp^na^erf c 2 " r |f%T r 

, p r r Zn(a-bb)-a.+r l j 

f»-tt) rcrfc iro t)'/* _P r 

4 K^)y, r| €Xp[-^(m.-»tt-r)][L<S rfc 

-ierdc^g^Jdr 

-^(ULH r 2e\pf^(m-rii+r)l\{erfr ZnC oM+^tc 

as-yosj J n^> 1 L D /J L 2(Dty^ 

- /> >-X, 2 >?fr-^t)-a 4 -r 7 J r 

6 r ^ c £fDt)'^ J ar 

- 3f£)- 3^r+^ 3 

+J (W'X^)[ ieric ^ +Cerfz ww i ] 

- 6(f)^e^-tV/c f^j + 
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The remaining integrals in the expression (E-18) for the fractional re- 
lease may be evaluated by integrating by parts, as follows: 



ol = 



du.= t~(^)exp[-$(m-hk>t-ij]d]r +expE^(™-r)bt'-rj]d)r ) 

<iv = er-fc — - Ar > 

^ { LnSy* 

and v-2(Dtpierfz . 



Several repetitions may be required before the pattern that develops 
becomes obvious. The result is 
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-l^: z . 



-Jr 



z(Dty 



+ LI /rv-W A-Okru -nrt- m)l ?/l iftL A To 2 nbY^S/xl ^ 2nCa~b£)± C 

T7v— /~r ld\ - yj^ z v v"? V/>CI 07 ^/ j l — 2?Dt)'te 



(E- 19) 
a.-r 



The similar integral, which would have dv = ierfc Q — ^ ± a — — dr, 

2 nTdF 

differs from Eq. (E-19) only in the power of the ierfc x term. If m 
were replaced by m + 1, the result [Eq„ (E-19)] could be used as it is, 
The remaining form of the integral to be evaluated by parts is 

(X~ y exp[- $-(na.-nkti-r)]> 
du.= -r(^)exp[p^(na.-M-hi'jjdt + exp [~jy( od-n&t+rjj dr > 

A y . er l A j 



and 



2hCo--kt)±& + r 

K- -(-v ' 



After several repetitions the result may be recognized as 




Upon applying limits to Eqs. (E-19) and (E~20), the resulting terms are 
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inserted into Eq. (E-18) to give 

f = 3(g)-3(#fH£T + 3($?' 

+ ^mj fa$] 

4- i **>& fen-lXa-ht) -al 

2(bk)^ 



+ 



-/■ 



/ y*h*£n fen+lIoL-kd-CL 1 

c Z(Dtyfe J 



+ 6 (^)Sf-0%-/)[k(^x^l /z ] 

m -z / U ., /i( Dt) '*■ 

+ t h te r-fc J 

- ( E -24) 

)9o “ 2. . 

+6(^V9|r / WJ^*p[-^^ ^rfc<&0§$±±- 



__ ',m / o r f r (Zn-lXd-bt)-<x 

L ert - ztkfi 



H" 



- (T ; e rfc 



CL 






- rer- 

1^7-3 



~-fc Z 

^-zft 



a-kt) - 0- 



+Q_ 

ZCDt)^ 



Otf 7 ^ 



-I2(§) *$ -l)7m-zl2(fXDtfj e,pl.f(nHXa-b§r € r^-^ m7z 

m=3 v y 

fe^ira-fcd-a- 7, 

(. en c 2rD0 )ft J’ 
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as the solution in full form for the fractional release of the tracer 
element from a spherical body. 

At t = 0, each term of Eq, (E-2 1) becomes zero; thus the frac- 
tional release is zero. At time t = a/b, which corresponds to com- 
plete evaporation of the solid, the fractional release should be 1, as 
follows : 




i -£c 2 . ^ ^ oj^dclJ/z 

= 3 - 3+1 = /• 



If, in Eq. (E-21), b is set equal to zero, which corresponds to no 
evaporation, then 



f 



.1 _ ^ X 0 - 7~ "I 

^.y ^ I , V- ^ 5 1 - (!X/ ,Z J ^ ° ~ (Crqj/J-j 



1 ^^ 03 / 1 } (lX)‘/z I j [L'-'n. J-oen^^/z J 



+ 
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+5(#J + £erfc&!& +ierf&i& + icriLAS*. 1 






( Dt)V: 



(W^ 



■(W-J 




- e(#)| [i : ^ _ i 

+ Q (#g [ftr/c (- - ; - ftrfc 




= ^(Wi^&r-fc. O + te/Vc^jv*] 



(E -23) 



-IZvtfL^&rfi O-cter/cT&n 



\"o> 




4 - ?(Dt) /z < YierfrJhib*- u. pr^ 424 ;, , rn-/)a. i 

T ^ ( <*J fci ^ - w*- 1 - l - n -W^ Ler ^rm7z i 



(m^J 



— f^(-^r}?. fiV^ — / k44 (n+Qck. _ ;**y<C(hrt)Q. j;^r(h-l)a' 



) 3 xoy-I^ o,'^, Y, /y. ,'3 , /* Cn~/)n 1 

~ Jb t) ,/s - ^ c c/ ' c (Dt)4; ^ d//"' C >rw'\/4 /• 




By letting n + 1 - m and n - 1 = p, one obtains 



Utir-fc 0 +i&Y-fc j^-/ 2 ~i~i~IZ(-^'l C%> tC~^~ 



(D tf* 



(-p-) z \l 3 er-fc 0 -i 3 er-f~ : 



+ 
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mo. 



<*V \^Z { ^fe,^ rT< - (cxy /z r £o tC (Dty*j 



_jC>/Dt\ y (*- er £r> P °~ -5 ;^ r 4 £23*. 1 
lf ^ \ 02- ) 1 ' ^ (DtY'2- ^-L. mrY^ i 



(og 



w=; 



W*J 



- 1 Q / D'C 'j^ /~ ^ I’^Prf'- — 0< »' ^(5f-Cr HQ = . . ”>7 pO- ~j 

If A r>n ) I ot tO Wro-yfe A-f-j O'-' 1 - (IX)‘'Z- ‘ ^ (ix) l/z J 



' ^ l IF/ tfe 



*W 



n=i 



3 (- 



Dtp r o : nr L^ n + *?rp»C. MQ- . o ? ,wZ. m,- jl<< :* r JZ _E°- 
oFj L^ cr ' ^ ^ ^ rt ^(Dt)"* ' p v, ter, '-(Dt)' / * 



..l \o/PL) f <? P /na_ _ ? /'VrA -££=.- — f ~ o 1 
< M a 2 -/ ^ u ' T{ -(Dty'z L rrc {ot)'z - - r '^ ° J 



P=' 

CL 



( E -24) 



lOfDtfzf^ ’ iZp^fr 



' ',^U/IL 

n-t 



t £ L 5 e r fz -j^r /z -f c 5 e rfc 0 - C 3 cr-fc 0 ] 



— 6 ("ga) [iGhfc 0 -h,2^z L-eri 



X, J2^ 1 

w*j 



l2(-^)L z erfc O ■ 



Because ierfc 






and 



.2 - 
l erfc 



0 = l/4, Eq. (E-24) becomes 



h Sff f k-‘ 



n=i 




( E -25) 



which may be recognized as Eq. (B-8). 
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E.3, Conversion to Dimensionless Form 



If the definitions (D-30) and (D-31) are now utilized, along with 
Eq. ( D-32), then Eq.(E-13) becomes 

/I 



'2&Z)i, exp ~ 



+er* 






^ Z exp{-np[hO-pzy r /o.]j fer-fc 

+ erfc2 na^t)-\+ r JZ\ 




Equation ( E -2 1) becomes 

-f = 3pT-3(^c)+(j3tf+- 5Z / H l ~/ 3 ' ,: X< :er fc^k + i^r-fc 2 ^k) 

- l2Z" z (^ erfc ^ -^ c2 0hr)- &Z(L z erfc^ 

- C z er-fc ) + 6 Zd-pzfi’er-fi: 4 ^ + c z erfc ) 

+3ZH^§f H (^'T^P 

+ i/vi-fc ] +- 



-9 3- 



)H 



( E -27) 



+3 3 [c r/c (2nt^-j3z)±i_ 

+ i»erfc &^^T . { \ 

4- 6 zSj-ir%-iX2”pT ,/z Y'exp [- 

-f ^<? rfc &±lXlzkl)zi 

ZT'/i J 

-f 6 X(l-Bz)Z(- l) m (2n0Z' / f%p[-n(3(n-iXi-/3z)][i^e r3c^2^0t± 

+ 6 Z(i-^(-07znp^r expt^iXi-pzM^er-Fc ^nz^tl 

m:z ' ‘ -L^er-fc &+'*£*>-> -f Z ~ 

+/2 Z 3/z £(-l) m (M ~^X^ n P >/P T e x pt^ n ~ l fi -frflJl^rXcfcz 1 ^ 1 ]^-- 

- 12 >/z T exp 

-rerfz^^p-], 



m-3 

and Eq. (E-25) becomes 



-f= G r' /2 [n-' ,& +2|!, Cer^c^J- 3Z • (e- 28) 



E.4. Approximate Solution for the Fractional Release 

Now, if the restrictions ^ <0.01 and (3t< 0.5 are made, then 
Eq. (E-27) may be closely approximated by 



f = 3(3C-3( f S'cf-f^z) 5 +3'r' / Y / -/ s '^ 'fcj^7 z - 
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i-"erfc 

Zj I C- 

(E -29) 

+6 Z^'irt^-^pZ^expj^ nfln-iXl-ptf] Zeri S 2 ^'/^ - 

n--Z ^ L 



-6Z(l^z)S(-ir(2n^TeypE^ Fer kfrfjte - 

m-z 222 




m=3 



of which the n = 1 term is most significant. 



Thus, f- 3p'Z-3(pzf+(pzf+5Z' /z (l-^)ierfc^p' /z ' 

- &Zi\rfc £p*+ 6 Z(l-p?)iW-kp^-lzf /z C 3 erfi 

43T' /2 ('-/3 z)lp m ^f\ppf' r er f c ^ 

i^rfc (E- 

- 6 Z(I-(SZ)Z(- r^n-^fZcrfc =§p* 

-12 z 3/z Z(-irh-zXHr'^TZ m er2c=M ,/z ■ 

hn =3 ^ ^ 
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Application of the inversion formulae for i m e rfc ( -x) - -developed in 

Appendix C--to Eq. (E-30), upon expansion to a sufficient number of 

/ a T 1/2 \x 

terms, reveals that terms in ^ ) cancel each other. There- 

upon the approximate expression for the fractional release will reduce 
to 

-f = 3pz-3((SzT+(f3'zf-3Z(i-f3z)-3Z 

+ 6 zfi-/Sz)ierfc -§f’' z + 12 r C( i-pz)i3er-fc 

+ Tl rlc j£' z 

(E -3 1) 

+6Z(i-pz)g(«n^T^er{cgf'* 

+ GzS(m + ^r(^Ti n e^c ^ 

m=i <• £ 

+ iz z zk 2(t»HX«T(gf*) m 

In general, the last three summations of Eq. (E-31) may be neglected 
in comparison with the remaining terms 

f= 3fz(l-pz)+(f2Z.f-3T-3Z(l-f3Z) 

3GZ k ( l-fz) ierfc-^%/ZZ(i-^) CZrfc 



(E -32) 
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NO ME NC E A T U RE 

a Characteristic dimension at initial conditions; half thickness for 
a slab; radius for a sphere (cm) 
b Rate of boundary movement (cm/sec) 

C Concentration of the tracer element (atoms/cm ) 

Cg Initial concentration of the tracer (atoms/cm ) 

D Diffusion coefficient for the tracer atoms (cm^/sec) 
f Fractional release of the tracer element 

Q Number of tracer atoms 

Qq Initial value of Q 

r Space coordinate in spherical geometry (cm) 
t Time from start of evaporation/diffusion (sec) 

x Space coordinate in slab geometry (cm) 

(3 Dimensionless e vaporation/diffusion coefficient [Eq. (II-8)] 

£ Integration variable (cm) 

T Dimensionless time [Eq. ( II - 7 ) ] 
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mission, nor any person acting on behalf of the Commission: 

A. Makes any warranty or representation, expressed or 
implied, with respect to the accuracy, completeness, 
or usefulness of the information contained in this 
report, or that the use of any information, appa- 
ratus, method, or process disclosed in this report 
may not infringe privately owned rights; or 

B. Assumes any liabilities with respect to the use of, 
or for damages resulting from the use of any infor- 
mation, apparatus, method, or process disclosed in 
this report. 

As used in the above, ’’person acting on behalf of the 
Commission” includes any employee or contractor of the Com- 
mission, or employee of such contractor, to the extent that 
such employee or contractor of the Commission, or employee 
of such contractor prepares, disseminates, or provides access 
to, any information pursuant to his employment or contract 
with the Commission, or his employment with such contractor. 
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